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Abstract. Siboglinids, previously referred to as pogono- 
phorans, have typically been divided into two groups, frenu- 
lates and vestimentilerans. Adults of these marine proto- 
stome worms lack a functional gut and harbor 
endosymbiotic bacteria. Frenulates usually live in deep. 
sedimented reducing environments, and vestimentiferans 
inhabit hydrothermal vents and sulfide-rich hydrocarbon 
seeps. Taxonomic literature has often treated frenulates and 
vestimentiferans as sister taxa. Sclerolinim has traditionally 
been thought to be a basal siboglinid that was originally 
regarded as a frenulate and later as a third lineage of 
siboglinids. Monilifera. Evidence from the 18S nuclear 
rDNA gene and the 16S mitochondrial rDNA gene pre- 
sented here shows that Sclerolinum is the sister clade to 
vestimentiferans although it lacks the characteristic mor- 
phology (ie., a vestimentum). The rDNA data confirm the 
contention that Sclerolbuun is different from frenulates, and 
further supports the idea that siboglinid evolution has been 
driven by a trend toward increased habitat specialization. 
The evidence now available indicates that vestimentiferans 
lack the molecular diversity expected of a group that has 
been argued to have Silurian or possibly Cambrian origins. 


Introduction 


Siboglinids were formerly called pogonophorans and in- 
clude two groups of marine protostomes, frenulates and 
vestimentiferans, that are commonly referred to as beard- 
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worms and tubeworms, respectively. Both groups lack a 
functional gut as adults and rely on endosymbiotic bacteria 
for nutrition. They have a closed circulatory system and 
possess a metamerized tail region called the opisthosoma. 
Vestimentiferans are distinguished from frenulates by the 
presence of a vestimentum, a winged region near the ante- 
rior of the organism. Both taxa occur in reducing environ- 
ments and typically are found at depths below several 
hundred meters. Due to the limited availability of samples 
and the difficulty of retrieving live specimens, several as- 
pects of their biology (e.g., reproduction, physiology) are 
still poorly understood. Vestimentiferans, in general. have 
been better studied than frenulates because they are key- 
stone species in eastern Pacific hydrothermal vent habitats 
and in Pacific and Caribbean seeps. 

The taxonomic literature concerning frenulate and vesti- 
mentiferan siboglinids has a colorful and confusing history. 
One taxonomic scheme recognizes frenulates (aka pogono- 
phorans sensu stricto) and vestimentiferans as distinct phyla 
(Jones, 1985). Alternatively, vestimentiferans have also 
been recognized as a class within the phylum Pogonophora 
(Jones, 1981: Ivanov, 1994). Others place frenulates and 
vestimentiferans within the phylum Annelida (Land and 
Nørrevang, 1977; Kojima et al., 1993; Bartolomaeus, 1995: 
McHugh, 1997: Rouse and Fauchald, 1997; also see South- 
ward, 1988). The latter hypethesis has been supported by 
recent morphological (Rouse and Fauchald, 1995, 1997), 
embryological (Young et al., 1996; Southward, 1999), and 
molecular analyses (Kojima et al., 1993: McHugh, 1997; 
Black et al., 1997: Kojima, 1998; Halanych et al., 1998). To 
further complicate matters, a ranked classification scheme 
has produced different names for the same clade of organ- 
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isms. Vestimentiferans have been called Vestimentifera 
(Jones, 1981), Obturata (Jones, 1981: Southward, 1988; 
Southward and Galkin, 1997), and Afrenulata (Webb, 
1969). Frenulates have been called Pogonophora (Jones, 
1985), Frenulata (Webb. 1969), Perviata (Southward, 
1988), and originally Siboglinidae (Caullery, 1914). 

Hereafter we apply the following nomenclature: (1) Ves- 
timentifera are equated with Obturata and Afrenulata; (2) 
Frenulata are equated with Perviata and Pogonophora 
(sensu Jones, 1985); (3) Monilifera is a third monogeneric 
clade that includes Sclerolinum,; and (4) Siboglinidae refers 
to the clade that includes Vestimentifera, Frenulata, and 
Monilifera. We recognize that the term “Pogonophora™ is 
more commonly used and that rules of priority tor nomen- 
clature do not apply to higher taxa. However. we have opted 
to use the term “Siboglinidae™ throughout this manuscript to 
emphasize that this group of organisms represents derived 
annelids (McHugh, 1997: Rouse and Fauchald. 1997). We 
restrict the term “pogonophoran™ to common usage. 

Even among siboglinids, there has been one group, 
Sclerolinum, that has been particularly problematic in terms 
of phylogenetic position. Unlike most frenulates that live in 
the mud, Scleroliman species can live on decaying organic 
material like wood or rope made from natural fibers (Webb, 
1964a: Southward, 1972). This taxon was originally con- 
sidered a member of the frenulate family Polybrachtidae 
(Southward, 1961), but Webb (1964b), mainly citing differ- 
ences in the postannular region, argued that Sclerolinum 
could not be ascribed to either of the two orders (Theca- 
nephria and Athecanephria) of siboglinids recognized at the 
time (vestimentiferans had not been discovered yet). He 
erected a new family. Sclerolinidae. that he states should 
“have order rank.” Ivanov (1991) more formally recognized 
the unique nature of Sclerolinum, and in 1994 he proposed 
that Frenulata (= Perviata), Monilifera (= Sclerolinidae), 
and the Vestimentifera be regarded as three taxa with equal 
rank (7.e., classes within the phylum Pogonophora). Addi- 
tionally, Ivanov (1994) further suggested that Monilifera are 
allied to the Vestimentifera on the basis of the common 
absence of several characters (e.g., spermatophores, teloso- 
mal diaphragm. metasoma preannular and postannular re- 
gions) relative to the Frenulata. Southward (1999) sug- 
gested that Monilifera might be similar to the ancestral 
siboglinid form, thus predicting that it should occupy a basal 
position in siboglinid phylogeny. Distinguishing between 
these hypotheses on the placement of Sclerolinum will 
allow us to test the notion of Black et al. (1997) that habitat 
preference or specificity may be an important factor in 
siboglinid evolution. If Black et al. are correct, Sclerolinum 
is expected to occupy a position between frenulates and 
vestimentiferans (which may be consistent with Ivanov’'s 
ideas), and not a position basal to the frenulate-vestimen- 
tiferan clade. 


To date, molecular studies that include siboglinids have 
either focused on vesumentiferans (Williams er al., 1993: 
Black et al., 1997: Kojima et al.. 1997: Halanych et al., 
1998) or have addressed siboglinid origins (Winnepen- 
ninckx et al.. 1995a; Kopma et al.. 1993; Kojima, 1998; 
McHugh, 1997). Most studies have included only one frenu- 
late representative. Although Black er al. (1997) included 
two “frenulate” siboglinids. one of these, the Loihi worm, 
was undescribed. Additionally. several 18S sequences were 
reported in a symposium contribution (Halanych ef al., 
1998) for which page limitations did not permit detailed 
analyses or explanation. Herein we extend these previous 
analyses by increasing the sampling of frenulates, including 
Sclerolinum, and using novel 18S rDNA and 16S rDNA 
data. The present findings support the notion that habitat 
requirements have been important in siboglinid evolution, 
Additionally, trenulates are sister to a Sclerolinum-vest- 
mentiferan clade, the latter of which showed limited diver- 
sity suggestive of a recent radiation within the clade. 


Materials and Methods 
Taxa employed 


Table | lists the species analyzed and GenBank accession 
numbers for the rDNA sequences used in this study. The 
frenulate and vestimentiferan operational taxonomic units 
(OTUs) included in this study represent all of the currently 
recognized genera available to the authors. The addition of 
closely related species within a genus would have increased 
OTUs without increasing the phylogenetic signal for the 
issues under examination and were therefore excluded. For 
example, there are no nucleotide differences observed in the 
18S rDNA of Escarpia spicata (Guaymas Basin) and E. 
laminata (Florida Escarpment). Limiting the number of 
OTUs also reduced computation time, allowing for more 
thorough analyses. Unless otherwise noted, collection local- 
ities correspond to those given in Black et al. (1997). 
Siboglinum ekmanii S. fiordicum, and — Sclerolinum 
brattstromi were collected near Bergen, Norway, and iden- 
tified by Eve Southward, Marine Biological Association of 
the United Kingdom. Identification of the frenulates Spiro- 
brachia and Polybrachia were made by Eve Southward on 
the basis of animal and tube morphology. Both specimens 
were collected by TVGrah from the Aleutian Trench 
(57°27.394'N. 148°00.013’W) at a depth of 4890 m on the 
German research vessel Sonne. 

The non-siboglinid annelid OTUs tor the 18S data were 
chosen to represent a diversity of lineages for which se- 
quences were available. The arthropod (Artemia) sequence 
was designated as the most distant outgroup for rooting 
purposes. Based on both morphology (e.g., Eernisse et al., 
1992) and molecular studies (e.g., Halanych ef al., 1995; 
Winnepenninckx er al., 19953; Aguinaldo et al., 1997; 
Eernisse, 1997), arthropods are clearly outside of the proto- 
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TABLE 1 
Taxa used in rDNA analyses 
GenBank Accession" GenBank Accession? 
Organism ISS rDNA 16S rDNA Organism 18S rDNA 16S rDNA 
Pogonophora Chaetopterida 
Frennlata Chaetopterus variopedatus U67324 
Galathealinum brachiosum AF168738 AF315040 Hirudinea 
Polybrachia sp. AF168739 AF315037 Haemopis sanguisuga X91401“ 
Siboglinum fiordicum GB X79876° Hirudo medicinalis AF315058 
Siboglinum fiordicum AF3 15060 AF315039 Oligochaete 
Siboglinum ekmani AF315062 AF315038 Enchytraeus sp. Z83750" 
Spirobrachia sp. AF168740 AF315036 Phyllodocida 
Vestimentitera Glycera americana U19519% 
Escarpia spicata AF168741 AF315041 Polynoidea 
Escarpiid n. sp. AF315053 Lepidonotopodium fimbriatum AF315056 
Lamellibrachia barhami AF168742 AF315043 Branchipolynoe symmytilida AF315055 
AF315044 Sabellida 
AF315045 Sabella pavonina U67144 
AF315047 Tubificidae 
Oasisia alvinae AF168743 AF315052 Tubifex sp. AF315057 
Rideeia piscesae AF168744 AF315048 Echinra 
AF315051 Ochetostoma ervitirogrammon X79875" 
AF315054 Urechis sp. AF315059 
Ridgeia piscesae GB X79877° Sipuncula 
Riftia pachyptila AF168745 AF315049 Phascolosoma granulatum X79874° 
AF315050 Nemertea 
Tevnia jerichonana AF168746 AF315042 Lineus sp. X79878" 
Monilifera Mollusc 
Sclerolinum brattstromi AF315061 AF315046 Scutopus ventrolineatus X91977! 
Annelida Priapulida 
Alvinellidae Priapulus caudatus X80234 
Paralvinella palmiformis AF168747 Arthropod 
Artemia salina X01723" 


“Unless otherwise noted, sequences were obtained in this study. 


P Sequence from Winnepenninckx et al. 


(199543). 


° Sequence from Nadot and Grant (unpublished). 


3 Sequence from Kim er al. (1996). 
° Sequence from Halanych et al. (1995). 


f Sequence from Winnepenninckx er al. (1996). 


© Sequence from Winnepenninckx er ai. 
h Sequence from Nelles et al. (1984). 


(1995): 


stome worm radiation. Because siboglinids are not closely 
related to molluscs and because of rate helerogeneity prob- 
lems within the Mollusca, only a single representative (the 
aplacophoran Scutopus) was used. Due to alignment limi- 
tations, oulgroups employed in the 16S analyses—a leech, 
an oligochaete, two polynoid polychaetes, and an echiu- 
rid—were more limited (see Table 1). Because different 
investigators collected the data at different times, there was 
not a ]:] correspondence in OTUs between data sets. We 
felt it more important to present all the relevant data rather 
than trim taxa from the data sets. The aligned data sets are 
available at the journal’s Supplementi’s page (http:// 
www.mbl.edu/BiologicalBulletin/VIDEO/BB.video.html) 
and at TREEBASE (http://phylogeny.harvard.edu/treebase). 


Data collection 


Total genomic DNA was extracted using a modified 
hexadecyl-trimethyl-ammonium bromide (CTAB) protocol 
(Doyle and Dickson, 1987). The entire 18S nuclear rDNA 
gene was amplified via PCR (polymerase chain reaction), 
using the universal metazoan oligonucleolide primers 1|8e 
and 18P (Halanych et al.. 1998). A region of the 16S 
mitochondrial rDNA was amplified using 16Sar-5’ and 
16Sbr-3’ primers (Palumbi, 1996). Each 50 ul reaction 
consisted of about 50 ng of template DNA, 0.5 uM of each 
primer, 2.5 mM MgCl, 200 uM dNTPs, 5 ul of manufac- 
turer's 10X reaction buffer, and 1.5 U Taq polymerase 
(Promega lne., Wisconsin). Cycling profiles were as fol- 
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lows: 18S—initial denaturation at 95 °C for 3 min. 35 
cycles of amplification (denaturation at 95 °C for | min, 
annealing at 50 °C for 2 min, extension at 72 °C for 2 min 
30 s). and a final extension at 72 °C for 5 min; 16S—ainitial 
denaturation at 94 °C for 2 min, 40 cycles of amplification 
(denaturation at 94 °C for 30 s, annealing at 46 °C for 30 s, 
extension at 72 °C for ] min), and a final extension at 72 °C 
for 7 min. PCR products were purified using the QIAEX II 
gel extraction kit (Qiagen Inc., California). Approximately 
60 ng of purified PCR product was used in sequencing 
reactions according to the manufacturer's instructions (FS 
Dye Termination Mix or Big Dye, Applied Biosystems Inc., 
California), The reaction profile was 25 repetitions of de- 
naturation at 94 °C for 30 s, annealing at 50 °C for 15 s, and 
extension at 64 °C for 4 min. Dye-labeled fragments were 
separated by electrophoresis on a Perkin Elmer ABI 373A 
or 377 DNA sequencer. Both strands of the PCR product 
were sequenced. In addition to the PCR primers, the oligo- 
nucleotide primers used for sequencing are given in 
Halanych er al. (1998) or Hillis and Dixon (1991). The 
sequences were assembled and verified using the AutoAs- 
sembler and Sequence Navigator programs (Applied Bio- 
systems Inc., California). The terminal primer regions were 
not included in the sequences submitted to GenBank or in 
the phylogenetic analyses. 


Phylogenetic analyses 


Sequence alignment was produced with a Clustal W 
program (Thompson et al., 1994) and subsequently cor- 
rected by hand using the protostome secondary structure 
models available through the Ribosomal Database project 
(http://rdp.cme.msu.edu/html/). Regions that could not be 
unambiguously aligned (e.g., divergent loop domains) were 
excluded from analyses. Tree reconstructions were imple- 
mented with the PAUP* 4.0b4b2 program (Swolford, 
2000), and MacClade 3.06 (Maddison and Maddison, 1992) 
was used for character and tree analyses. Neighbor-joining 
(NJ), parsimony, and maximum likelihood (ML) analyses 
were performed and yielded similar results. In the interest of 
brevity, results and discussion will focus on ML analyses. 

NJ trees were reconstructed under Jukes-Cantor, Kimura- 
2-parameter, Tamura-Nei, general-time-reversible, and log/ 
det models. All except log/det were examined under equal 
rates of among-site rate variation using the empirically 
derived gamma shape parameter. a, of 0.3 (see Swofford er 
al., 1996, tor summary of different assumptions used in 
these models). A Kishino-Hasegawa (1989) likelihood eval- 
uation of the resulting topologies revealed no significant 
differences between models for either the 16S or the 18S 
data. Kishino-Hasegawa evaluations estimated a six-substi- 
tution-type rate matrix for which nucleotide base frequen- 
cles were set to empirical values and œ was estimated. NJ 
bootstraps consisted of a log/det correction (model was 


arbitrarily chosen) with 1000 iterations. Parsimony analyses 
consisted of heuristic searches with 100 random sequence 
additions and tree-bisection-reconnection (TBR) branch 
swapping. Transitions (Ti) and transversions (Tv) were 
given equal weighting. ML evaluation of parsimony topol- 
ogies was the same as for NJ topologies. One thousand 
iterations were used for parsimony bootstrap analyses. 
When using likelihood to search for the “best” tree (as 
opposed to evaluating given trees). computation time was 
limiting. Therefore, we used a nucleotide model with two 
substitution types where the Ti/Tv ratio was set to the value 
estimated for the best parsimony tree (empirical base fre- 
quencies were used). ML searches were heuristic with 10 
random sequence-addition replicates. ML bootstraps em- 
ployed the “Faststep” option with 100 iterations. 


Results 


The 18S rDNA data set consisted of 26 OTUs and 1935 
nucleotide positions. Of the 1614 nucleotide positions that 
could be unambiguously aligned, 34.6% (559 positions) 
were variable and 18.7% (303 positions) were parsimony 
informative. Figure 1 shows the single best likelihood tree 
(Ln likelihood = —8260.55148) recovered. All search 
methods in all analyses found a monophyletic siboglinid 
clade (bootstrap support was 298% for all methods). Res- 
olution within the vestimentiferan clade, as well as between 
annelid groups, was poor, however. The moniliferan 
Sclerolinun brattstromi falls out with the vestimentiferan 
taxa in all analyses (bootstrap = 98%). The remaining 
frenulates form a distinct sister-clade to the Sclerolinum- 
vestimentiferan clade with =99% bootstrap support. 

Resolution among annelid taxa and within the vestimen- 
tiferans was poor due to the lack of phylogenetic signal. 
Because this paper does not focus on the annelid radiation, 
we did not try to enhance resolution among all annelid taxa. 
However, we did attempt to boost the signal within the 
vestimentiferan clade by employing a less inclusive taxo- 
nomic alignment. For metazoan 18S sequences, inclusion of 
broader taxonomic diversity can often create larger regions 
of ambiguous alignment that should not be included in 
analyses, due to poor assumptions about positional homol- 
ogy. Thus by reducing the taxonomic breadth examined, the 
phylogenetic signal can potentially be increased by a “bet- 
ter” alignment (Halanych, 1998). Unfortunately, even when 
just the siboglinids were aligned, little genetic diversity was 
observed, and the vestimentiferan taxa were still poorly 
resolved (not shown). The exception was Lamellibrachia 
barhami, which was consistently placed as the most basal 
vestimentiferan. Table 2 shows the logdet/paralinear dis- 
tances (below diagonal) and absolute distances (above di- 
agonal) for this less-inclusive, siboglinid-only alignment (in 
which most divergent domains could be unambiguously 
aligned). Even though the distance values for the siboglinid- 
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Figure 1. Results of 18S rDNA phylogenetic anatyses. The single best likelihood tree (Ln likelihood = 


—§$260.55148) found. Analysis details are given in the text. Maximum likelihood bootstrap values of =50% are 
given in bold. Parsimony (italic) and neighbor joining (underlined) values are also given for the major nodes of 
interest (values for other nodes were omitted in the interest of space). Branch lengths are drawn proportional to 
the inferred amount of change along the branch (scale shown). 


only alignment are only slightly greater than the full align- 
ment values, the greatest distance within vestimentiferans 
was only 0.02 (with a maximum of 25 nucleotide differ- 
ences), revealing that there was very little 18S genetic 
diversity within this group. 

The 16S rDNA data set consisted of 24 OTUs, each with 
497 nucleotide positions. Of the 465 nucleotide positions 
that could be unambiguously aligned, 60.4% (281 positions) 
were variable and 47.7% (222 positions) were parsimony 
informative. The reconstructed topology (Ln likelihood = 
~ 3967.21062), Figure 2, was qualitatively similar to the 
18S topology. Siboglinids are divided into two major lin- 


eages: vestimentiferans plus the moniliferan Sclerolinum 
brattstromi (bootstrap support 83% for ML and 100% for 
NJ and parsimony) and a frenulate sister-clade (bootstrap 
support 294% in all analyses). Again, S. brattstromi was 
basal to the vestimentiferans. In a departure from the 18S 
analyses, Riftia pachyptila, not Lamellibrachia barhami, 
often fell out as the most basal vestimentiferan. However, 
this was never supported by >54% bootstrap support; ML 
analyses that excluded the non-siboglinid outgroups re- 
vealed that the base of the Vestimentifera was poorly re- 
solved with 16S data. A comparison of genetic divergence 
values (Table 3) indicates that there was limited genetic 
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TABLE 2 


Pairwise distances for the siboglinid-only 18S rDNA data set; absolute distances above diagonal and log/det distances below diagonal 








l 2 3 4 5 6 7 8 9 10 ll 12 13 14 

l Spirebrachia — 109 113 87 ie 13] 124 [22 125 124 131 is t2] 120 

2 Polybrachia 0.07 -— 2 104 138 137 Ibe 140 140 140 147 142 140 136 

3 Galathealinian 0.07 0.01 _ 106 142 141 142 143 143 143 150 145 143 139 

4 Siboglinum ekmani 0.05 0.06 0.06 — 116 117 113 110 113 110 121 ji 107 pie 

5 Siboglinum fiordicum 0.08 0.09 0.09 0.07 — 5 140 136 136 140 143 138 134 l. 

6 Siboglinum fiordicum GB 0.08 0.08 0.09 0.07 0.00 -— 143 l. 139 143 146 141 [37 140 

7 Escarpia 0.08 0.09 0.09 0.07 0.09 0.09 -- 7 i4 10 l 6 13 3i 

8 Ridgeia 0.08 0.09 0.09 0.07 0.09 0.09 0.00 -- 8 7 Ly 4 ]2 a2 

9 Ridgeia GB 0.08 0.09 0.09 0.07 0.09 0.09 0.01 0.00 — 14 21 ll 19 38 

10 Oasisia 0.08 0.09 0.09 0.07 0.09 0.09 0.0] 0.00 0.01 — 20 J 14 32 
ll Riftia 0.08 0.09 0.09 0.07 0.09 0.09 0.01 0.01 0.01 0.01 — 16 29 39 
]2 Teviiia 0.08 0.09 0.09 0.07 0.09 0.09 0.00 0.00 0.01 0.00 0.01 —— ll 30 
13 Lamellibrachia 0.07 0.09 0.09 0.07 0.08 0.09 0.01 0.01 0.01 0.01 0.01 0.01 — 25 
14 Sclerolinum 0.08 0.09 0.09 0.07 0.09 0.09 0.02 0.02 0.02 0.02 0.02 0.02 0.02 — 


variation within vestimentiferans (£0.11 log/det distance; a validated by the present data. Moreover, we found that 


maximum of 47 nucleotide differences). Sclerolinum brattstromi is closely allied to the vestimenti- 

As for the frenulate clade, neither 18S or 16S supported a ferans, and does not occupy a position basal to a frenulate- 
monophyletic Siboglimun; bul because only two Siboglinum vestimentiferan clade, confirming Ivanov’s (1991; 1994; 
species were examined, additional taxa are needed to verify the Ivanov and Selivanova, 1992) ideas that moniliferans oc- 


status of this frenulate taxon. Additionally, we performed cupy a position intermediate between vestimentiferans and 
Kishino-Hasegawa (1989) likelihood evaluation for both genes frenulates. 

to test the monophyly of the frenulate and vestimentiferan- Southward (1993) also suggested a possible evolutionary 
Sclerolinum clades. To this end. we used the constraints option link between Sclerolinwn and vestimentiferans. This con- 
in PAUP* 4.0b4b2 to conduct parsimony heuristic searches tention is confirmed by the present analysis, as well as a 
(specifics same as above) to find the best trees that were recent morphological cladistic analysis (Rouse, 2001). Us- 
consistent and inconsistent with the monophyly of these ing 44 morphological characters coded for all recognized 
clades. Both the 16S and the 18S data significantly support the — sjboglinid genera, Rouse found support for the monophyly 
monophyly of both groups (18S trenulates—average ML score — of Frenulata, Vestimentifera. and the Scleroliniun-vestimen- 
supporting monophyly = —8244.69. non-monophyly score = tiferan clade. However, our use of nomenclature differs 
~$278.135. P value < 0.01; 16S trenulates—monophyly = from Rouse with regard to the term Monilifera, which he 
— 3894.889. non-monophyly = —3927.49, P value < 0.005: applies 10 the Sclerolinuim-vestmentiferan clade. Because 


ISS vestimenhferan-Sclerofinium—monophyly = — 8244.69. this term was originally (Ivanov and Selivanova, 1992) 
non-monophyly = —827].922. P value < 0.05; 16S vestimen- applied to only Sclerolinum, and because of the morpho- 
tiferan-Sclerolinuuan—monophyly = —3894.889, non-mono- 


logical differences from vestimentiferans, Rouse’s use of 
the term will inject confusion into the literature. Although 
we acknowledge that Monilifera, as defined here, is redun- 


phyly = —3911.802. P value < 0.05). 


Discussion dant with the generic name Sclerolinuim, several aspects of 

The monophyly of siboglinids (aka, Pogonophora sensu siboglinid evolution and taxonomy are in need of additional 

lat) is supported by morphological (Southward. 1988, study. Thus. we have chosen not to name this clade until 
1993: Rouse and Fauchald, 1995; Rouse, 2001). embryo- more is understood about siboglinid evolution. 

logical (Southward. 1999), and molecular (Winnepenninckx The placement of Sclerolinum was especially interesting 


et al., 1995a: Black et al., 1997; McHugh, 1997: Halanych in the context of the evolution of habitat preference. Previ- 
et al., 1998, this study) evidence. Thus. in agreement with ous studies of vestimentiferans (Black et al., 1997), clams 
others (Southward, 1988, 1999; Ivanov, 1994; McHugh, (Peek et al.. 1997). mussels (Craddock et al.. 1995), and 


1997). we see no support for the recognition of vestimen- shrimp (Shank er al., 1999) reveal that vent-endemic organ- 
titerans and frenulates as having fundamentally different isms are related to. and possibly derived from, species 


body plans (.e., “phyla” sensu Jones, 1985). The assertion associated with hydrocarbon seeps that occur near subduc- 
made by Webb (1964b) and later by Ivanov (1991, 1994) tion zones and continental margins. Furthermore, recent 
that Sclerolinuwm was notably different from frenulates is observations (Feldman et al., 1998: Baco et al., 1999; Distel 
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Resnlts of 16S rDNA phylogenetic analyses. The best likelihood tree (Ln likelihood = 


—3967.21062) found. Another tree with a Ln likelihood score of —3967.25739 was fonnd in the same search. 
The trees differed in relationships within the Ridgeia clade. Analysis details are given in the text. Maximum 
likelihood (ML) bootstrap values of =50% are given in bold. Parsimony (italic) and neighbor joining (under- 
lines) valnes are also given for the major nodes of interest (values for other nodes were omitted in the interest 
of space). In the ML bootstrap analysis, Lamellibrachia and Sclerolinwn formed a clade in 55% of the iterations. 
That is not shown above becanse it is incompatible with the “best” ML tree. Branch lengths are drawn 
proportional to the inferred amount of change along the branch (scale shown). 


et al., 2000) reveal that several symbiont-bearing clams. 
vestimentiferan tubeworms, and mussels can survive on 
rotting organic material, such as wood or a whale carcass. 
The moniliferan S. brattstromi and related species (e.g., S. 
javanicum, S. minor, and S. major) are typically found 
growing on decaying organic material such as wood or rope 
(Webb, 1964a, b; Southward, 1972: Ivanov and Selivanova, 
1992). Other members of the genus, (e.g., S. sibogae and S. 
magdalenae) lived buried in mud (Southward, 1972). These 


habitat preferences suggest that affinity for a mud or silt 
habitat was ancestral in siboghinids, allowing us to speculate 
that a pattern of evolution from low-oxygen, sedimented 
habitats to decaying organic material to hydrocarbon seeps 
to hydrothermal vents has occurred within the Sclerolinum- 
vestimentiferan clade. 

Although neither the 18S nor the 16S data clearly resolve 
relationships within the Vestimentifera, the cytochrome c 
oxidase subunit I (COI) data of Black et al. (1997) show 
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TABLE 4 


Percent of significant tests when comparing relative substitution rates 
between the two major siboglinid clades 


Number Significant Percenl 
Tem type* of tests results significant 
Between trenulates and 
veslimenliferan- 
Sclerolinum clade 
{SS rDNA 480 331 69.0 
t6S rDNA 350 47 13.4 
Within frenulates 
tSS rDNA 150 53 55.3 
16S rDNA 50 l 2.0 
Within vestimentiferan- 
Scleroliniuim clade 
ISS rDNA 280 SO 28.6 


16S rDNA 455 iS 2.9 


Results of relative rates tests basedon an HKY plus gamma model in the 
HyPhy program. The program is distributed by S. Muse, Department of 
Slalistics, North Carolina State University. 

* The 18S comparisons employed all Lophotrochozoan outgroups. 


seep tubeworms to be basal to vent tubeworms (but see 
Wiliams et al, 1993). This pattern in the evolution of 
habitat preference roughly proceeds from less reducing to 
more reducing (greater sulfide and methane availability) 
environments. A similar evolutionary trend was observed in 
bathymodiolid mussels (Craddock et al., 1995; Distel et al., 
2000). Examination of additional taxa is needed to verify 
whether this is a general trend in the evolution of vent and 
seep taxa. 

All molecular studies to date (Williams er al., 1993: 
Black et al., 1997; and Tables 2 and 3) reveal that vesti- 
mentiferans exhibit very limited molecular diversity for a 
group suggested to be several hundred million years old. 
This lack of diversity may be due to a slowdown in the rate 
of molecular change (7.e., nucleotide substitution) in vesti- 
mentiferans, a recent common origin for extant vestimen- 
tiferans, or possibly both. For the present 18S rDNA se- 
quences, vestimentiferans appear to have expertenced a 
significant molecular slowdown relative to the frenulates or 
other protostome taxa (Table 4: as judged using an Hk Y 
plus gamma correction model in the HyPhy software pack- 
age distributed by S. Muse, Department of Statistics, North 
Carolina State University). With 16S data, only 13.4% of 
tests between frenulates and members of the vestimenti- 
feran-Sclerolinum clade were significant. Although this 
Value is not statistically significant, it is a greater percentage 
than is found within either group (~3%), suggesting that a 
limited rate discrepancy may exist. Similar rate disparities 
Were not observed for COI data (Black et al. 1997), but 
only one frenulate was used in the comparison. Nonetheless, 
we concluded that present-day vestimentiferans constitute a 
young evolutionary group. 


In contrast, previous interpretation of Silurian tubeworm 
fossils (Little er al., 1997) as vestimentiferans suggested 
that these worms constitute an ancient animal lineage. It is 
possible that the Silurian tubeworm fossils represent an 
earlier offshoot from an ancient siboglinid lineage, but this 
will be impossible to test as the fossils lack the necessary 
soft-tissue preservation. Additionally, we note that many 
wormlike invertebrates make tubes. For example. some 
alvinellid polychaetes observed during our recent expedi- 
tion to vents along the Southern Eastern Pacific Rise (32°S., 
]00°W) occupted tubes with diameters comparable to the 
tubes of mature Aiftia pachyptila. Many of the alvinellid 
tubes were partially overgrown by sulfide chimneys, and 
thus were effectively “fossilized.” Although we are not 
convinced of the interpretation of Silurian fossils as repre- 
sentative of an extant lineage of vestimentiferans, we should 
point out that specimens from the Cretaceous are convinc- 
ing (Little er al, 1999). In contrast, all the hydrothermal! 
vent-endemic taxa that have been examined with appropri- 
ate molecular tools appear to be from relatively recent 
radiations (i.e. <100 MY; Black er al.. 1997: Peek et al., 
1997; Shank et al., 1999; McArthur and Koop, 1999; but see 
McArthur and Tunnicliffe, 1998. for possible exceptions). 
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